Long-lived spin entanglement induced by a spatially correlated thermal bath 
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We investigate how two spatially separated qubits coupled to a common heat bath can be entan- 
gled by purely dissipative dynamics. We identify a new dynamical timescale associated with the 
lifetime of the dissipatively-generated entanglement and show that it can be much longer than either 
the typical single-qubit decoherence time, or the timescale on which a direct exchange interaction 
can entangle the qubits. We give an approximate analytic expression for the long-time evolution of 
the qubit concurrence and propose an ion trap scheme in which such dynamics should be observable. 
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I. INTRODUCTION 

Entanglement is the hallmark of correlations in quan- 
tum theory and has come to be seen as a precious re- 
source essential for many quantum information process- 
ing protocols An entangled state can show corre- 
lations stronger than those allowed classically, however, 
these correlations are often extremely fragile. Interac- 
tions with the environment surrounding any quantum 
system tend to cause a rapid loss of quantum coherence 
2] , generally leading to entanglement within the system 
being destroyed on a short timescale 0. Various tech- 
niques have thus been developed to protect entangled 
states from their surroundings, such as constructing deco- 
herence free subspaces [3, [1[ , dynamical decoupling 0, 0] , 
and exploiting the quantum Zeno effect [8J. 

As well as avoiding decoherence, entanglement must 
also be generated. This can be achieved by a number of 
means; for example, by harnessing intrinsic system cou- 
plings @ , through projective measurements [l(| , or via a 
quantum bus [ll| . Recently, it was shown that entangle- 
ment between a pair of two-level systems (qubits) can in 
fact be generated by the same processes that are usually 
considered to be detrimental, if the qubits are allowed 
to interact with a common bath [l^, EH . This offers a 
potential way to explore the interplay between coherent 
and incoherent multi-qubit dynamics with a significantly 
reduced level of external system control. 

In general, immersing a pair of otherwise non- 
interacting spins in a common heat bath will give rise 
to two terms in the subsequent master equation: a uni- 
tary, Hamiltonian-like Lamb-shift term 0, [HJ leading to 
coherent spin evolution, and a dissipative term [2J, Il3l . 
both of which may generate spin entanglement [12], Il3l . 
[Til . [TEl [iH ] . While it is known that in the idealised case of 
unseparated qubits entanglement induced by the dissipa- 
tive term can persist indefinitely [l7| . surprisingly little 
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attention has been paid to the dynamics of its generation 
and decay in the more realistic setting of finite inter-qubit 
separation. In this experimentally relevant case, it is not 
clear on what timescale the generated entanglement per- 
sists, or even whether its observation is feasible at all. 
It has been shown that the entangling capability of the 
Lamb-shift is highly sensitive to the inter-qubit separa- 
tion [l4| , with dissipative processes destroying any gener- 
ated entanglement more rapidly as separation increases. 
Furthermore, it has recently been shown that the level of 
entanglement generated between two harmonic oscilla- 
tors suffers from a similar critical dependence on oscilla- 
tor separation Hence, a comparison of the timescales 
associated with conventional decoherence dynamics to 
those for dissipatively-induced spin-entanglement gener- 
ation and decay is needed. 

In this article, we address the above issues by studying 
the dynamics of bath-induced entanglement in the con- 
text of the two-spin-boson model, which has wide appli- 
cations in the solid state and elsewhere [l9l |20| . We show 
that for small but finite spin separation, the timescale on 
which dissipativcly-induccd spin entanglement survives 
can be far larger than the corresponding single spin de- 
coherence time. It also exceeds the timescale on which 
entanglement induced by either a direct exchange inter- 
action or the Lamb-shift persists. In particular, we ob- 
tain an approximate analytic expression for the long-time 
dynamics of the two-spin concurrence and from this de- 
termine its survival time. We suggest an ion trap reali- 
sation of our model and demonstrate that observation of 
the generated entanglement should indeed be feasible. 



II. MASTER EQUATION 



We consider two spatially separated, identical, non- 
interacting spin qubits, each subject to a static field of 
strength A/2 in the x-direction, and coupled to a com- 
mon bath of harmonic oscillators. The Hamiltonian is 
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given by 



H = 
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with ft = 1. Here, cr™ is a Pauli operator acting on the 
nth qubit (n — 1, 2; i = a;, y, z), _B„ is the bath operator 
coupling qubit n to the bath, while LUk is the angular fre- 
quency, and b, (bk) the creation (annihilation) operator 
of the bath mode of wave vector k. 

To investigate the dynamics of the reduced two-spin 
density operator p we assume that the qubit-bath cou- 
pling is weak compared to A and follow the stan- 
dard Born-Markov and rotating-wave approximation ap- 
proach 0. This relies on a perturbative expansion in 
the system-bath coupling strength, and the assumption 
that the bath instantly re-thermalises after any interac- 
tion with the qubits. Setting (Hi) s = 0, we obtain a 
Schrodinger-picture master equation of the usual form 



d P (t) 
dt 



-i[H s + H L3 ,p(t)]+V(p(t)), 



(2) 



valid to second order in Hi. Here, the Lamb-shift pro- 
vides a Hamiltonian-like contribution and is of the form 



H LS = A(a\ + a 2 x ) + B(a l z a 2 z + a 1 ^), 



(3) 



where we have omitted a constant. The precise expres- 
sions for A and B depend on the details of the qubit-bath 
interactions and will be given in section Hill For now we 
will comment on their likely qualitative effects. The first 
term of H^s simply renormalises the static field strength 
due to the presence of the bath modes. The entangling 
capability of the second term deserves some attention 
since it represents an induced interaction between the 
two qubits, as has been explored in detail in Ref. [3]. 
However, we will show that in general this entanglement 
decays at a rate far quicker than that generated by the 
dissipator, here given by 



£>(/>(*)) =E lnM(A m (L0)p(t)Al(L0) 
n,m— 1 

-{4(c)A m (c),p(i)}), (4) 



uj n,m=l 

~ 2 



with frequency summation over the eigenvalue differ- 
ences of Hs (u> — ±A) and corresponding eigenoperators 
A„(±A) = (1/2)(^ T ^). 



III. BATH CORRELATION FUNCTIONS 

The key quantities in the present discussion are the 
Fourier transforms of the bath correlation functions 
7mnM = dse luJS (Bl(s)B m ), for which we need a 
specific form for B n . As usual in the spin-boson model, 



we consider linear coupling between the qubits and the 
coordinate of each bath mode [2(| such that 
B n(s) = E k (9k b l eluJkS +9k* b ke- tWkS ), with coupling con- 
stants g%. Note that with this form of B n our assump- 
tion of (Hi)b = leading to Eq. ([2]) is justified. An 
important aspect of this work is that the qubits have an 
explicit spatial separation and thus g\ ^ g\. To make 
this evident, we consider our spins to be separated by a 
distance d along the z-axis such that g\ 



1 g^idkcosO /2 



and g\ = g ke - ldkcos() / 2 j where 9 is the polar angle mea- 
sured against the z-axis in fc-space, and \g\ \ — \gV\ = gk- 
Taking the bath to be in thermal equilibrium, we find 



712(0; 



dse^^| fffc | 2 (A( Wfe )e^ s e 

-°° k 

+ {N(ujk) + i) e -^ s e - ifcdcose ), 



ikd cos 6 



(5) 



where N(ujk) = [exp^fe/feeTs) — is the thermal oc- 
cupation of mode k, Ub is Boltzmann's constant, Tb the 
temperature of the bath, 712(0;) = 721(0;) (once the sum- 
mation is performed), and 7nn(o;) is obtained by setting 
d = 0. 

Defining the bath spectral density to be J(u>) = 
Sfc lfffc| 2 <5(o; — u>k) pjl, and an inverse dispersion relation 
(assumed isotropic) k = k(ui), we take the continuum 
limit of the summation over k above to find 



712(0;) = /(k(o;)cZ)7ii(o;), 



(6) 



with 711 (+ A) = 722(+A) = (N(A) + 1)70 and 
7n(-A) = 722 (- A) = A(A) 7o , where 7o = 2ttJ(A) 
is the single-spin decoherence rate at zero temperature. 
Here, f(x) describes the bath's spatial correlations and is 
determined by its dimensionality D(= 1, 2, 3). For D = 1 
we have f(x) = cos(a;); for D = 2, f(x) = Jq(x), where 
Jo is a Bessel function of the first kind; and for D = 3, 
f(x) = sinc(a;). We thus write 



7 i 2 (±A) = (l-<5)7n(±A), 



(7) 



where 1 — 8 captures the degree of correlation between 
the baths seen by each qubit, becoming unity at d = 
<5 = (completely correlated) and, for D > 1, zero 
as d — > 00 (5 — > 1, completely independent). When 
D = 1, the Markovian assumption constrains 712 to be 
periodic with respect to d. However, we will concen- 
trate here on the limit where d is small enough such that 
S « (n(A)d) 2 /2D, for all D. 

The strength of the two terms in the Lamb-shift Hamil- 
tonian are given by combinations of the Hilbert trans- 
forms of the bath correlation functions. They are found 
to be 

A = 2 J o J(u)coth(u/(2k B T B )) ( A2 ^ 2 ) d^ (8) 



and 



B= I J(u))f(n(u)d) ( A2 ) doj m 
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where principal values are assumed. Note that for 
system-bath coupling in 2 or 3 dimensions B — > as 
the qubit separation is increased to infinity, expressing 
the fact that an uncorrelated bath can not give rise to 
any coherent coupling between the qubits. In contrast, 
A contains no distance dependence since it represents a 
renormalisation of the single-qubit energy levels in each 
spin, independent of any bath correlations. 

The relative strength of the coherent terms in the evo- 
lution, A and B, compared to the strength of the dis- 
sipative terms, given by 7n and 712, dictates whether 
the bath is capable of generating entanglement through 
the Lamb-shift [l4[ ■ Evaluation of the relevant integrals 
involved necessitates a specific form of the bath spectral 
density. In this article we shall focus on the entanglement 
generated through dissipative processes and as such we 
may leave A and B unevaluated. In fact, we shall show 
in Seciton lVll that the dissipatively induced entanglement 
can persist for times far larger than entanglement gener- 
ated through the Lamb-shift, regardless of its strength. 

IV. STATE DYNAMICS 

Rather than working directly with the reduced density 
operator p, it is more instructive to work in terms of a 
16 dimensional vector a, which is a generalisation of the 
Bloch vector for two-qubit states. It is constructed by 
flattening the matrix whose elements ay satisfy 

1 3 

p=-2_,a ij a\®a 2 j , (10) 

i,i=0 

where <7q = <Tq = /. The traceless property of the Pauli 
matrices ensures that = (aja'j), and conservation of 
probability demands that «oo = 1. To describe the evo- 
lution of our system, we consider the eigensystem of the 
Liouvillian super-operator C defined (unconventionally) 
by a(t) — Ca{t) 1 where the linearity of the transforma- 
tion between p and a ensures that the dynamics gener- 
ated by C is equivalent to that of Eq. @. Clearly, a state 
a. initially in an eigenstate of £, say on (with eigenvalue 
A;), will evolve according to a.(t) = one Xlt . Hence, a 
general state evolves such that 

15 

a(t)=Y,aiaie Xlt , (11) 

1=0 

where the coefficients ai are determined from the initial 
conditions, and the sum runs over all eigenstates of C. 

We are primarily interested in the long-time dynamics 
of our system. Hence, it makes sense to search for eigen- 
values of C with small (or zero) real parts since, assuming 
these parts are all negative, the corresponding eigenvec- 
tors will contribute towards the total state Eq. pTjl on 
the largest timescale. We evaluate the full eigensystem of 
C analytically, though this leads to combersome expres- 
sions which we shall not give here. However, we can iden- 
tify a number of important features for our subsequent 



analysis. For any non-zero qubit separation (<5 ^ 0) there 
exists a single eigenvector, oiq, with zero eigenvalue, and 
all other eigenvalues have negative real parts. We there- 
fore associate c*o with the thermal state since it is that to 
which all states tend towards as t — > 00. Of the remaining 
15, there is a single eigenvalue, Ai (with corresponding 
eigenvector ot\), which vanishes as 5 — > at all temper- 
atures. Expanding the exact expression to first order in 
5 and second order in -/V(A) gives the simple form 

Ai = -(1 + 3JV(A))*7 . (12) 

It is also possible to show graphically that Ai varies 
approximately linearly with N(A) at all temperatures. 
For reasons that should become clear, we shall refer to 
the eigenvector corresponding to the null eigenvalue, ao, 
and the eigenvector corresponding to the vanishing eigen- 
value, ai, as our eigenvectors of interest. The eigenvec- 
tor corresponding to the thermal state, a , is expressible 
solely in terms of R = (l + 2iV(A))" 1 = tanh(A/2fc B T.B), 
and is given by 

a : {a 00 , a i, an, 0^22} = {1, R, R 2 , 0}, (13) 

where the notation implies that the eigenvector has the 
elements specified, and that 022 = 033, am = a w, with 
all other elements being zero. From a numerical analysis 
of ol\ we find that for 5 <C 1, to a very good approxima- 
tion ol\ can also be expressed just in terms of R, with 
corrections being of the order of 6: 

oti : {a 00 , a i, an, a 2 2} ~ {0, R, 1 + R 2 , 1}, (14) 

where the notation is the same as in Eq. (|13[) . 

There are two further notable eigenvalues, A2 and A3 = 
A2. Once again, expanding the exact eigenvalues to first 
order in <5, and to first order in R about R = 1, gives the 
expression 

A 2 = A3 = -(1/2)70(1 -R+ 2S- SR) - iA. (15) 

Note that these eigenvalues have vanishing real parts only 
in the limit that both the separation and temperature go 
to zero (S — ► and R — ► 1). In either the zero tempera- 
ture limit (R — > 1) or the zero separation limit (<5 — * 0) 
the corresponding eigenvectors, 0L2 and 0:3 = ctj, are 
given by 

£* 2 : {a 2, a 3, ai 2 , ai 3 } = {i, 1, i, 1}, (16) 

where once again the notation implies the eigenvector has 
the elements specified but this time with «20 = — ^02, 
a?3o = — «03, o^i = — «12 and 0:31 = — «i3, and all other 
elements being zero. 

This leaves us with 12 eigenvalues to consider. Plotting 
their real parts as a function of R it becomes clear that 
they are all ~ — 70/-R, for all values of S. As such, the 
corresponding eigenvectors contribute towards the total 
state significantly only for times t < ■fQ 1 regardless of the 
temperature. These eigenvalues and eigenvectors shall be 
referred to as those with / > 3. 



4 



With the relative size of the real parts of the various 
eigenvalues in mind, we see that for small enough 6, at 
times sufficiently greater than 7^ , a general state can 
be approximated by 



a(t) = a Q + aiaie Xlt + a 2 (a 2 e X2t ± ct^e"^') 



(17) 



where we have normalised olq, set clq = 1, and set 
03 = ±a 2 to ensure positivity of the corresponding den- 
sity operator. The coefficient a\ is found, by setting 
t = in Eq. (11]), to be a\ = (A — i? 2 )/(3 + R 2 ), where 
A = (cr 1 • <t 2 ) = an + a 22 + CK33- For separable pure 
states A represents a scalar product of single qubit Bloch 
vectors. Positivity of the corresponding density operator 
limits the range of A to —3 < A < +1. 

To gain insight into the general features of a state de- 
scribed by Eq. (fT?)) . it is useful to write down the corre- 
sponding density operator in the zero temperature and 
separation limit, in which there are no dccoherent pro- 
cesses due to the real parts of the relevant eigenvalues 
vanishing. Using Eq. (fTT))) we find 

pit) =(1 + ai) IHUHL - oi |*"> X <*"L 

+ a 2 V2(e- lAt \*-) w m\ x ± e lAt ITT> a (*-| B ) , 

(18) 

where = (1/V2)(\U) X UOJ = (l/V2)(|Ti), - 

UT) J is the Bell singlet, and ITT), = (1/2)(|TT} 2 + IU} Z + 
|TI) Z + |lf) z )- The subscripts x and z refer to the rel- 
evant basis, with our qubits defined with respect to z. 
Written in this way, we can see that our state is a coher- 
ent mixture of the singlet and the state |TT) X - At zero 
temperature, the state \\\) x is the ground state since 
it minimises the energy associated with the the system 
Hamiltonian and energy can not be absorbed from the 
environment. Also, at zero separation, the Bell singlet is 
stable since it is the totally anti-symmetric state, while 
the Hamiltonian is totally symmetric [22]. Therefore, in 
a combination of these limits coherent mixtures of these 
states are stable. However, the states are at different 
energies which gives rise to the exponential factors in 
Eq. HID. 



V. ENTANGLEMENT DYNAMICS 

We are now in a position to say that, provided our 
qubits are sufficiently close together such that 5 <C 1, for 
times t > Jq 1 contributions from eigenvectors with / > 3 
and their associated dynamics will have all but vanished, 
and to a good approximation (and a better approxima- 
tion as time increases), our two-qubit state will be well 
described by Eq. (fT?) . To quantify the resulting entangle- 
ment dynamics we use Wootters concurrence [23j . which 
ranges from for separable states to 1 for maximally 
entangled states. It can be shown numerically that the 
concurrence of a state described by Eq. (fT?) depends only 




FIG. 1: (Color Online) Concurrence of the initially separable 
state |1" J.) in either x, y, or z (A = —1) as a function of time 
(scaled by Ai) and R = tanh(A/2fcsTB) calculated from the 
full Liouvillian, though ignoring the Lamb-shift. 



very weakly on the magnitude of a 2 , and in the zero tem- 
perature and separation limit has completely vanishing 
dependence. We may therefore set a 2 = to derive a 
simple concurrence expression, and from Eqs. (|13p . (|14|) 
and (JTTJ) find 



[R 2 - 1){R 2 + 3) + (R 2 - A) (3 - R 2 )e 



Ait 



° = maX [ 2WT3^ ' ° 

valid (with increasing accuracy) for timescales t > 7 . 
The legitimacy of setting a 2 = will be demonstrated 
towards the end of this section, where comparisons with 
numerics using the full eigensystem of C are made. 

Whether the bath is capable of inducing spin entangle- 
ment, and if so how it subsequently decays, is now clear. 
Firstly, no entanglement is generated unless 



A < 



5R 2 



R 2 



(20) 



in agreement with Ref. [171 ] . Secondly, provided this in- 
equality is satisfied, we see from Eq. ([T9]) that the induced 
entanglement will decay exponentially until (R 2 — A)(3 — 
R 2 )e Xlt = -(R 2 -l)(R 2 + 3) is satisfied, after which time 
the entanglement is zero. This occurs at 



tc = [ATI ln 



{R 2 -A)(R 2 -3) 



{R 2 + 3)(i? 2 -1) 



(21) 



as demonstrated in Fig. [TJ Note that Eqs. (fT5[) and 
(|21[) are valid for a range of temperatures, however, in 
view of the inequality of Eq. (f2l3) , we will focus on small 
temperatures since these maximise the amount and life- 
time of any induced entanglement. Interestingly, in the 
limit of vanishing temperature (R — > 1), the entangle- 
ment reaches zero only asymptotically. Also, as the qubit 
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separation d — > the level of entanglement becomes a 
function of the initial state onl y a nd t c — > oo, i.e. the 
steady-state becomes entangled [17| . In general, t c varies 
as (S"fo)~ 1 for given R and A, hence it can be length- 
ened by increasing the ratio A/Tg, or by decreasing the 
separation d. 

We would naturally like to know which initially separa- 
ble states result in the largest generated concurrence. For 
fixed temperature and spin separation the only parame- 
ter to consider is A, and from inspection of Eq. (| 19(1 we 
see that it should be minimised. This corresponds to an 
initial state that is as anti-symmetric as possible; hence, 
A is minimised by anti-aligning the single spin Bloch vec- 
tors, giving A = — 1 for pure states. Such a state corre- 
sponds to |f |) in x, y or z. As the states become more 
mixed, the Bloch vectors decrease in length and A — * 0. 
Interestingly, even a maximally mixed state (A = 0) can 
become entangled provided A > 2 cothT 1 (y/ 5/ 3) IcbTb- 
In general, we expect an initially separable state to reach 
its maximum entanglement after a time ~ 7^ 1 , typical 
of the decay of eigenvectors a; for I > 3. Note that if 
we do not restrict our initial state to a separable state, 
A is minimised by the Bell singlet, for which A = —3. As 
we have mentioned, in the limit that the qubit separa- 
tion goes to zero the singlet is able to maintain its full 
entanglement for all times. 

To illustrate these points, in the main part of Fig. © 
we plot the time evolution of the concurrence for vari- 
ous initial states, calculated both from Eq. (fl"9|) (dashed 
lines) and numerically using the full Liouvillian (solid 
lines), though here ignoring the Lamb-shift. As claimed, 
on timescales > Jq 1 ~ 5 in the scaled time units, the 
entanglement dynamics is well approximated by the an- 
alytic form. Note also that neglecting the eigenvectors 
a.2 and 0:3 has had no discernible effect on the accu- 
racy of Eq. (fT!)]) on the timescales it is expected to be 
valid. For the initially separable states (A = —1,0), we 
clearly see that the entanglement reaches its maximum 
on a timescale > Jq 1 ~ S, decaying subsequently at a 
rate ~ Sjo. For the Bell singlet (A = —3), the analytic 
approximation becomes almost exact since this state is 
simply a linear combination of the vectors cxq and a.\ in 
Eqs. (fl"3|) and (]14f) . Also shown is the behaviour of the 
Bell state (l/-s/§)(|Ti>, + UT)„), far which A = 1. Un- 
like the singlet, all other Bell states have the maximum 
possible value of A, and as such lose their entanglement 
rapidly. 

We could consider the maximally mixed state (A = 0) 
as being the infinite temperature thermal state since it 
represents a state for which thermal fluctuations com- 
pletely overcome any external fields. We see that as 
the bath "cools" this state towards the thermal state 
at Tbi it does so via entangled states. Of course, af- 
ter a time t c the state of the qubits becomes separable 
once more, and will eventually reach the thermal state 
at Tb- There is in fact a well defined condition for Tb 
and (initially prepared) qubit temperature Tq such that 
the bath has the ability to entangle the qubits. In the 
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FIG. 2: (Color Online) Main: Concurrence as a function 
of (scaled) time calculated analytically (dashed lines, valid 
for t > 7(7 1 ) and numerically (solid lines). We consider 
four initial states, the Bell singlet (l/v^XIU) - lit}) (blue, 
A = —3), the pure state |fj.) (red, A = —1), the max- 
imally mixed state p = (1/4)/ (green, A = 0) and the 
Bell state (1/V2)(|f j.) + |lf)) (black, A = 1). Parameters: 
S = 0.05, R = 0.9. Inset: Behaviour of the initial state |f|) 
where, in the numerical calculations, the Lamb-shift and ex- 
change interactions have been included at an arbitrarily cho- 
sen strength B = £ = l/(2|Ai|). 



limit of small bath and qubit temperature this condition 
becomes 9b — 0q > ^m3, where 8b — (A/2/s;,Tb) and 
9 Q = (A/2k b T Q ). 



VI. ENTANGLEMENT GENERATED 
THROUGH THE LAMB-SHIFT 

It is important to consider the role of the Lamb-shift 
Hamiltonian Hl$, thus far ignored. Within the limits of 
our derivation, namely that A is large and hence rota- 
tions about x are so rapid that the y- and z- directions 
are effectively equivalent, its form is determined by sym- 
metry. With this form it commutes with our eigenvectors 
of most interest, a and ol\. It can also be shown that 
the eigenvectors a.2 and a.3 are eigenoperators of Hls, 
with eigenvalues 2(A + B) and —2(A + B), respectively. 
Therefore, the Lamb-shift Hamiltonian can influence only 
the eigenvectors with I > 3, the imaginary parts of their 
eigenvalues, and the imaginary parts of A2 and A3. Hence, 
despite the fact that H^s can entangle our spins, its effect 
is restricted to a timescale of order 7^ (after which the 
other eigenvectors have decayed) regardless of its ampli- 
tude, and it will therefore have no effect on the long-time 
entanglement dynamics or the analytic expressions we 
have derived. 

In fact, we can further account for the effect of a direct 
spin exchange interaction simply by adding a term of the 
form He = ^cr 1 - a 2 into Hg in Eq. @, provided that the 
evolution it generates occurs on timescales much slower 
than the bath correlation time. This procedure is valid 
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in the regime of A 3> £, such that A sets the relevant 
frequency scale for the system-bath interaction. In this 
case, exactly the same conclusions hold as for the Lamb- 
shift term since olq and aj again commute with this form 
of interaction, and ol-i and 0:3 are also eigenoperators of 
He- Its influence will thus similarly be restricted to short 
timescales ~ ■ This point is illustrated in the inset of 
Fig. @ where we plot the analytically and numerically 
calculated concurrence of the initial state ||J.) in x, y or 
z, including both H^g and He with arbitrarily chosen 
strengths. As expected, their impact is seen only on a 
timescale ~ 7q~ , much shorter than that over which the 
dissipativcly induced entanglement survives. 

VII. EXPERIMENTAL REALISATION 



we find a revival time ~ 507^ . Assuming the wave- 
length associated with A to be approximately N / (A/uJt) 
in units of the ion spacing, we choose to address two 
neighbouring ions to give 8 = (A/ui t N) 2 /2D w 0.03. 

Furthermore, the temperature of the bath must be low 
enough such that the inequality of Eq. (j2TJ)) is satisfied, 
allowing a finite level of entanglement to be generated. 
By requiring R — 1/2 we find that the bath must have 
a temperature in the mK range for typical trapping fre- 
quencies of MHz. We conclude that after a timescale 
~ 7^ = (257raw t ) _1 a concurrence of C ~ 0.15 should 
be generated from the initial state |||) in x, y or z. It 
will subsequently decay by at least a factor of e _1 before 
the dynamics associated with the finite size of the bath 
becomes significant. 



While spin-boson models apply commonly in solids, 
more controlled realisations in ion traps have recently 
been proposed [24j| . We extend to the two-spin-boson 
model as follows: we consider a linear ion trap with the 
internal levels of a single ion representing each spin, cou- 
pled to the collective motion of N ions providing a (finite) 
bosonic bath. The two-spin-boson Hamiltonian is gener- 
ated by simultaneously addressing two ions with travel- 
ing waves which, in a linear ion trap, gives rise to an 
Ohmic spectral density J(uj) = (a/2)ui. The static field 
A is set by the laser-ion Rabi frequency. The strength 
of the system-bath interaction can be adjusted by vary- 
ing various experimental parameters such as the laser 
wavelength and ion mass. We work in the weak coupling 
regime and therefore set a — 0.1 Addressing the 

ions in the way described also induces an effective Ising 
interaction between the two ions due to the polaron rep- 
resentation 24j , which nevertheless disappears for a finite 
lattice with Ohmic spectral density. 

Owing to the finite size of the bath, the system evolves 
as if it were coupled to a continuum only for short times, 
after which a quantum revival is seen. For example, for 
N = 100 this revival occurs at approximately 2n/uJt, 
where oj t is the trapping frequency 24} . Hence, to ob- 
serve both the generation and subsequent decay of the 
dissipatively-induced entanglement we describe, this pe- 
riod must be larger than both ^q 1 and (c>7o) -1 . Since 
both 70 and 5 depend on A, this corresponds to a careful 
choice of A. With A = 25u>t (set by the Rabi frequency), 



VIII. SUMMARY 

We have shown that the mechanisms normally as- 
sociated with dissipative processes can lead to long- 
lived entanglement in non-interacting, spatially separated 
two-qubit systems. We have highlighted two important 
timescales. The first, shorter timescale 7 c j~ 1 is that with 
which we expect a single qubit to dephase. When a 
second qubit is introduced close to the first, we find 
dissipatively-induced entanglement is generated on this 
short timescale, and further that there is a second larger 
timescale (c^o) -1 on which the induced entanglement de- 
cays. Importantly, the influence of both the bath- induced 
Lamb-shift or a direct spin exchange interaction is still 
restricted to the original shorter timescale. Hence, the 
presence of a second qubit within the bath induces co- 
herences in the overall system state that can persist on 
timescales far larger than either the corresponding single 
qubit decoherence time, or timescales associated with the 
influence of direct exchange or the Lamb-shift. 



A. Acknowledgements 

We thank H. T. Ng for useful discussions. SB thanks 
the Royal Society and Wolfson Foundation. DPSM, AN 
and SB are supported by the EPSRC. 



[1] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
2000). 

[2] H.-P. Breuer and F. Petruccione, The Theory of Open 
Quantum Systems (Oxford University Press, 2002). 

[3] T. Yu and J. H. Eberly, Phys. Rev. B 66, 193306 (2002). 

[4] D. A. Lidar, I. L. Chuang, and K. B. Whaley, Phys. Rev. 
Lett. 81, 2594 (1998). 

[5] P. Zanardi and M. Rasetti, Phys. Rev. Lett. 79, 3306 



(1997). 

[6] G. Gordon and G. Kurizki, Phys. Rev. Lett. 97, 110503 
(2006). 

[7] F. F. Fanchini and R. d. J. Napolitano, Phys. Rev. A 76, 
062306 (2007). 

[8] S. Maniscalco, F. Francica, R. L. ZafRno, N. L. Gullo, 
and F. Plastina, Phys. Rev. Lett. 100, 090503 (2008). 

[9] B. Kraus and J. I. Cirac, Phys. Rev. A. 63, 062309 
(2001). 



7 



[10] S. D. Barrett and P. Kok, Phys. Rev. A 71, 060310(R) 
(2005). 

[11] A. Blais, R. Huang, A. Wallraff, S. M. Girvin, and R. J. 

Schoelkopf, Phys. Rev. A 69, 062320 (2004). 
[12] D. Braun, Phys. Rev. Lett. 89, 277901 (2002). 
[13] F. Benatti, R. Floreanini, and M. Piani, Phys. Rev. Lett. 

91, 070402 (2003). 
[14] D. Solenov, D. Tolkunov, and V. Privman, Phys. Rev. B 

75, 035134 (2007). 
[15] S. W. J. An and H. Luo, Physica A 382, 753 (2007). 
[16] T. Choi and H. J. Lee, Phys. Rev. A. 76, 012308 (2007). 
[17] F. Benatti and R. Floreanini, J. Opt. B 7, S429 (2005). 



[18] T. Zell, F. Queisser, and R. Klesse, Phys. Rev. Lett. p. 
160501 (2009). 

[19] A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. Fisher, 
A. Garg, and W. Zwerger, Rev. Mod. Phys. 59, 1 (1987). 

[20] M. Dube and P. Stamp, Int. J. Mod. Phys. B 12, 1191 
(1998). 

[21] G. D. Mahan, Many-Particle Physics (Plenum, 1990). 
[22] P. Zanardi, Phys. Rev. A 57, 3276 (1998). 
[23] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 
[24] D. Porras, F. Marquardt, J. von Delft, and J. I. Cirac, 
Phys. Rev. A 78, 010101(R) (2008). 



